Loss of phylogenetic and functional originalities of woody cerrado species in simulated extinction scenarios

Originality measures how different a given species is from all other co‐occurring species regarding either their phylogenetic history or functional traits. Since it is important to preserve the various aspects of diversity and original species carry more phylogenetic or functional information, originality may be used to assign conservation priorities. Our goal was to evaluate the relationships between phylogenetic and functional originalities, and their simulated losses under extinction scenarios based on abundance, fire tolerance and habitat preference. We placed 100 plots in a cerrado reserve located in central Brazil, sampled all woody plants species within the plots, measured 14 functional traits and measured fire history. We assembled a phylogenetic tree and a functional dendrogram, with which we calculated the originalities. Phylogenetic‐ and functional‐based originalities were correlated. However, the loss of functional originality was different from random extinctions on the abundance and fire tolerance scenarios, whereas the loss of phylogenetic originality was not. When compared with phylogenetic originality, functional originality brought more information to be used in conservation strategies because it was sensitive to differences in species abundance and fire tolerance. Thus, the extinction of rare or fire‐sensitive species would result in important functional changes due to loss of distinctive traits.     

Predicting chromatin organization using histone marks

Genome-wide mapping of three dimensional chromatin organization is an important yet technically challenging task. To aid experimental effort and to understand the determinants of long-range chromatin interactions, we have developed a computational model integrating Hi-C and histone mark ChIP-seq data to predict two important features of chromatin organization: chromatin interaction hubs and topologically associated domain (TAD) boundaries. Our model accurately and robustly predicts these features across datasets and cell types. Cell-type specific histone mark information is required for prediction of chromatin interaction hubs but not for TAD boundaries. Our predictions provide a useful guide for the exploration of chromatin organization.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0740-z) contains supplementary material, which is available to authorized users.     

Comprehensive miRNA sequence analysis reveals survival differences in diffuse large B-cell lymphoma patients

BackgroundDiffuse large B-cell lymphoma (DLBCL) is an aggressive disease, with 30% to 40% of patients failing to be cured with available primary therapy. microRNAs (miRNAs) are RNA molecules that attenuate expression of their mRNA targets. To characterize the DLBCL miRNome, we sequenced miRNAs from 92 DLBCL and 15 benign centroblast fresh frozen samples and from 140 DLBCL formalin-fixed, paraffin-embedded tissue samples for validation.ResultsWe identify known and candidate novel miRNAs, 25 of which are associated with survival independently of cell-of-origin and International Prognostic Index scores, which are established indicators of outcome. Of these 25 miRNAs, six miRNAs are significantly associated with survival in our validation cohort. Abundant expression of miR-28-5p, miR-214-5p, miR-339-3p, and miR-5586-5p is associated with superior outcome, while abundant expression of miR-324-5p and NOVELM00203M is associated with inferior outcome. Comparison of DLBCL miRNA-seq expression profiles with those from other cancer types identifies miRNAs that were more abundant in B-cell contexts. Unsupervised clustering of miRNAs identifies two clusters of patients that have distinct differences in their outcomes. Our integrative miRNA and mRNA expression analyses reveal that miRNAs increased in abundance in DLBCL appear to regulate the expression of genes involved in metabolism, cell cycle, and protein modification. Additionally, these miRNAs, including one candidate novel miRNA, miR-10393-3p, appear to target chromatin modification genes that are frequent targets of somatic mutation in non-Hodgkin lymphomas.ConclusionsOur comprehensive sequence analysis of the DLBCL miRNome identifies candidate novel miRNAs and miRNAs associated with survival, reinforces results from previous mutational analyses, and reveals regulatory networks of significance for lymphomagenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0568-y) contains supplementary material, which is available to authorized users.     

Uridine supplementation exerts anti-inflammatory and anti-fibrotic effects in an animal model of pulmonary fibrosis

Rationale

Pulmonary fibrosis is a progressive disease with only few treatment options available at the moment. Recently, the nucleoside uridine has been shown to exert anti-inflammatory effects in different animal models, e.g. in acute lung injury or bronchial asthma.

Method

Therefore, we investigated the influence of uridine supplementation on inflammation and fibrosis in the classical bleomycin model. Male C57BL/6 mice received an intratracheal injection of bleomycin on day 0 and were treated intraperitoneally with uridine or vehicle. The degree of inflammation and fibrosis was assessed at different time points.

Results

Uridine administration resulted in attenuated inflammation, as demonstrated by reduced leukocytes and pro-inflammatory cytokines in the broncho-alveolar lavage (BAL) fluid. Furthermore, collagen deposition in the lung interstitium was also reduced by uridine supplementation. Similar results were obtained in a model in which animals received repeated intraperitoneal bleomycin injections. In addition uridine inhibited collagen and TGF-ß synthesis by primary lung fibroblasts, the release of pro-inflammatory cytokines by human lung epithelial cells, as well as the production of reactive oxygen species by human neutrophils.

Conclusion

In summary, we were able to show that uridine has potent anti-inflammatory and anti-fibrotic properties. As uridine supplementation has been shown to be well tolerated and safe in humans, this might be a new therapeutic approach for the treatment of fibrotic lung diseases.     

Distinguishing low frequency mutations from RT-PCR and sequence errors in viral deep sequencing data

Background

RNA viruses have high mutation rates and exist within their hosts as large, complex and heterogeneous populations, comprising a spectrum of related but non-identical genome sequences. Next generation sequencing is revolutionising the study of viral populations by enabling the ultra deep sequencing of their genomes, and the subsequent identification of the full spectrum of variants within the population. Identification of low frequency variants is important for our understanding of mutational dynamics, disease progression, immune pressure, and for the detection of drug resistant or pathogenic mutations. However, the current challenge is to accurately model the errors in the sequence data and distinguish real viral variants, particularly those that exist at low frequency, from errors introduced during sequencing and sample processing, which can both be substantial.

Results

We have created a novel set of laboratory control samples that are derived from a plasmid containing a full-length viral genome with extremely limited diversity in the starting population. One sample was sequenced without PCR amplification whilst the other samples were subjected to increasing amounts of RT and PCR amplification prior to ultra-deep sequencing. This enabled the level of error introduced by the RT and PCR processes to be assessed and minimum frequency thresholds to be set for true viral variant identification. We developed a genome-scale computational model of the sample processing and NGS calling process to gain a detailed understanding of the errors at each step, which predicted that RT and PCR errors are more likely to occur at some genomic sites than others. The model can also be used to investigate whether the number of observed mutations at a given site of interest is greater than would be expected from processing errors alone in any NGS data set. After providing basic sample processing information and the site’s coverage and quality scores, the model utilises the fitted RT-PCR error distributions to simulate the number of mutations that would be observed from processing errors alone.

Conclusions

These data sets and models provide an effective means of separating true viral mutations from those erroneously introduced during sample processing and sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1456-x) contains supplementary material, which is available to authorized users.     

Identification and characterization of expressed retrotransposons in the genome of the Paracoccidioides species complex

Background

Species from the Paracoccidioides complex are thermally dimorphic fungi and the causative agents of paracoccidioidomycosis, a deep fungal infection that is the most prevalent systemic mycosis in Latin America and represents the most important cause of death in immunocompetent individuals with systemic mycosis in Brazil. We previously described the identification of eight new families of DNA transposons in Paracoccidioides genomes. In this work, we aimed to identify potentially active retrotransposons in Paracoccidioides genomes.

Results

We identified five different retrotransposon families (four LTR-like and one LINE-like element) in the genomes of three Paracoccidioides isolates. Retrotransposons were present in all of the genomes analyzed. P. brasiliensis and P. lutzii species harbored the same retrotransposon lineages but differed in their copy numbers. In the Pb01, Pb03 and Pb18 genomes, the number of LTR retrotransposons was higher than the number of LINE-like elements, and the LINE-like element RtPc5 was transcribed in Paracoccidioides lutzii (Pb01) but could not be detected in P. brasiliensis (Pb03 and Pb18) by semi-quantitative RT-PCR.

Conclusion

Five new potentially active retrotransposons have been identified in the genomic assemblies of the Paracoccidioides species complex using a combined computational and experimental approach. The distribution across the two known species, P. brasiliensis and P. lutzii, and phylogenetics analysis indicate that these elements could have been acquired before speciation occurred. The presence of active retrotransposons in the genome may have implications regarding the evolution and genetic diversification of the Paracoccidioides genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1564-7) contains supplementary material, which is available to authorized users.     

Identification of NF-κB and PLCL2 as new susceptibility genes and highlights on a potential role of IRF8 through interferon signature modulation in systemic sclerosis

IntroductionSystemic sclerosis (SSc) and primary biliary cirrhosis (PBC) are rare polygenic autoimmune diseases (AIDs) characterized by fibroblast dysfunction. Furthermore, both diseases share some genetic bases with other AIDs, as evidenced by autoimmune gene pleiotropism. The present study was undertaken to investigate whether single-nucleotide polymorphisms (SNPs) identified by a large genome-wide association study (GWAS) in PBC might contribute to SSc susceptibility.MethodsSixteen PBC susceptibility SNPs were genotyped in a total of 1,616 patients with SSc and 3,621 healthy controls from two European populations (France and Italy).ResultsWe observed an association between PLCL2 rs1372072 (odds ratio (OR) = 1.22, 95% confidence interval (CI) 1.12 to 1.33, P_adj = 7.22 × 10⁻⁵), nuclear factor-kappa-B (NF-κB) rs7665090 (OR = 1.15, 95% CI 1.06 to 1.25, P_adj = 0.01), and IRF8 rs11117432 (OR = 0.75, 95% CI 0.67 to 0.86, P_adj = 2.49 × 10⁻⁴) with SSc susceptibility. Furthermore, phenotype stratification showed an association between rs1372072 and rs11117432 with the limited cutaneous subgroup (lcSSc) (P_adj = 4.45 × 10⁻⁴ and P_adj = 0.001), whereas rs7665090 was associated with the diffuse cutaneous subtype (dcSSc) (P_adj = 0.003). Genotype-mRNA expression correlation analysis revealed that the IRF8 protective allele was associated with increased interferon-gamma (IFN-γ) expression (P = 0.03) in patients with SSc but decreased type I IFN (IFIT1) expression in patients and controls (P = 0.02). In addition, we found an epistatic interaction between NF-κB and IRF8 (OR = 0.56, 95% CI 0.00 to 0.74, P = 4 × 10⁻⁴) which in turn revealed that the IRF8 protective effect is dependent on the presence of the NF-κB susceptibility allele.ConclusionsAn analysis of pleiotropic genes identified two new susceptibility genes for SSc (NF-κB and PLCL2) and confirmed the IRF8 locus. Furthermore, the IRF8 variant influenced the IFN signature, and we found an interaction between IRF8 and NF-κB gene variants that might play a role in SSc susceptibility.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0572-y) contains supplementary material, which is available to authorized users.     

Postoperative radiotherapy in prostate cancer: Analysis of prognostic factors in a series of 282 patients

Aim

To assess the outcomes of patients treated with postoperative RT in relation to the possible prognostic factors.

Background

Postoperative radiotherapy (RT) has been proved to reduce the risk of biochemical recurrence in high-risk prostate cancer patients. Baseline prostate specific antigen (PSA), pathological Gleason score (GS), positive surgical margins, nodal status and seminal vesicle invasion are independent predictors of biochemical relapse.

Materials and methods

The clinical records of 282 patients who underwent postoperative RT were retrospectively reviewed. The prognostic value of postoperative PSA, preoperative risk class, nodal status, pathological GS, margins status, and administration of hormonal therapy (HT) was analyzed.

Results

Postoperative RT was delivered with a median dose to the prostatic fossa of 66 Gy (range 50–72) in 1.8–2 Gy/fraction. Median follow-up was 23.1 months (range 6–119). Five-year actuarial biochemical disease-free survival (bDFS) and overall survival rates were 76% and 95%, respectively. Higher bDFS was found for patients with postoperative PSA <0.02 ng/ml (p = 0.03), low preoperative risk class (p = 0.01), pN0 (p = 0.003), GS 4–6 (p = 0.0006), no androgen deprivation therapy (p = 0.02), and irrespective of surgical margin status (p = 0.10). Multivariate analysis showed that postoperative PSA and Gleason score had a significant impact on bDFS (p = 0.039 and p = 0.05, respectively).

Conclusions

Postoperative RT with a dose of 66 Gy offers an acceptable toxicity and an optimal disease control after radical prostatectomy in patients with different risk features. A postoperative PSA >0.02 ng/ml could be considered as a prognostic factor and a tool to select patients at risk for progression.     

Cloning,expression, and physicochemical characterization of a new diheme cytochrome <Emphasis Type="Italic">c</Emphasis> from <Emphasis Type="Italic">Shewanella baltica</Emphasis> OS155

The 16-kDa diheme cytochrome c from the bacterium Shewanella baltica OS155 (Sb-DHC) was cloned and expressed in Escherichia coli and investigated through UV–vis, magnetic circular dichroism, and ¹H NMR spectroscopies and protein voltammetry. The model structure was obtained by means of comparative modeling using the X-ray structure of Rhodobacter sphaeroides diheme cytochrome c (Rs-DHC) (with a 37% pairwise sequence identity) as a template. Sb-DHC folds into two distinct domains, each containing one heme center with a bis-His axial ligation. Both secondary and tertiary structures of the N-terminal domain resemble those of class I cytochrome c, displaying three α-helices and a compact overall folding. The C-terminal domain is less helical than the corresponding domain of Rs-DHC. The two heme groups are bridged by Tyr26 in correspondence with the shortest edge-to-edge distance, a feature which would facilitate fast internal electron transfer. The electronic properties of the two prosthetic centers are equivalent and sensitive to two acid–base equilibria with pK _a values of approximately 2.4 and 5, likely corresponding to protonation and detachment of the axial His ligands from the heme iron and a pH-linked conformational change of the protein, respectively. Reduction potentials of −0.144 and −0.257 V (vs. the standard hydrogen electrode), were determined for the C- and N-terminal heme groups, respectively. An approach based on the extended Debye–Hückel equation was applied for the first time to a two-centered metalloprotein and was found to reproduce successfully the ionic strength dependence of E°′.